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ABSTRACT 

The  work  of  Lewis  and  Walecka  on  the  structure  of  the  giant  reso- 
nance in  Carbon  12  and  Oxygen  16  has   seemed  to  indicate  that  the 
single -particle,    shell-model  theory  gives  better  results  than  collective 
models   --in  particular  the  Goldhaber  -  Teller  and  Steinwedel  -  Jensen 
models.      The  transverse  form  factor  as  a  function  of  momentum 
transfer  is  observed  experimentally  to  fall  to  a  minimum  and  rise 
again  in  the  range  from  zero  to  0.  6  f        .      This  dip  is  predicted  by  the 
Brown  particle-hole  theory  model  but  is  not  predicted  by  either  of  these 
collective  models.      The  purpose  of  this  paper  is  to  develop  a  collective 
model  which  could  possibly  give  this  dip  in  the  form  factor.      Further 
work  is  necessary  to  determine  if  this  model  can  actually  do  so.     A 
semiclassical  treatment  is  used  to  find  the  matrix  elements  for  the 
assumed  model.      The  transverse  form  factor  is  then  calculated  for  two 
nuclear  charge  density  functions   --  constant  density  and  the  Fermi 
distribution. 
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1.       Introduction. 

Electron  scattering  is  a  very  useful  tool  for  the  study  of  nuclear 
structure,    since  the  cross-section  for  electron  scattering  from  a 
nucleus  can  be  related  to  the  basic  geometric  details  of  the  structure 
of  the  nucleus.      The  incoming  electron  is  assumed  to  interact  with  the 
nucleus  by  means  of  the  electromagnetic  field.      This  interaction  is 
known,    whereas  the  nuclear  interaction  (strong  force)  is  not. 

The  parameters  of  electron  scattering  are  the  initial  momentum   k 
and  initial  energy  E  of  the  incoming  electron,    which  are  controlled  by 
the  experimenter;  and  the  final  momentum   k     and  the  final  energy  E    , 
which  are  measured.      The  cross -section  can  be  found  as  a  function  of 
the  momentum  transfer    q      ;   k    -    k     from  the  electron  to  the  nucleus 
and  of  the  excitation  energy,    which  can  be  found  from  E,    E    ,    the 
scattering  angle,    and  considerations  of  the  nuclear  recoil.      One  can  then 
vary   q   for  a  fixed  excitation  energy  (i.  e.  ,    for  an  excitation  of  the 
nucleus  to  a  particular  level)  ,     and  get  the  cross-section  as  a  function 
of   q   (where    q  =   I  q  I  )  for  that  level.      This  cross- section  is  related  to 
the  Fourier  transform  of  the  charge  and  current  densities  of  the  nucleus. 
Thus,    we  can  get  information  about  the  charge  and  current  densities  of 
the  nucleus  from  electron  scattering. 

A  so-called  "dipole  state"  is  a  feature  of  the  electron  scattering 
cross-section  for  most  nuclei.     The  dipole  state  is  a  state  strongly 


excited  by  photons  of  the  appropriate  energy  (usually  10  -  25  Mev).      It 
shows  up  as  the  "giant  dipole  resonance",    which  is  a  resonance  peak 
in  the  cross -section  as  a  function  of  energy.      The  dipole  state  is   so 
named  because  it  has  a  large  electric  dipole  matrix  element  with  the 
ground  state.     As  a  result  of  this  last  fact,    these  states  lead  to  strong 
photon  absorption  cross -sections  which  are  orders  of  magnitude  larger 
than  those  of  other  multipoles. 

If  nuclear  recoil  is  ignored,    the  cross-section  for  electron 
scattering  can  be  written  in  the  general  form  [2]  : 


da  8ria2      k' 


duj  4  k 

q 


2kk' 

+  — 
q 


—      2kk'    cos2    (9/ 2)    f    2    (q) 

q 

sin2   (9/ 2)     [  q2  +  2kk '   cos2    (9/2)  ]    i*  (q) 


where  q    is  the  four -momentum  transfer  and  9    is  the  angle  between  the 

incoming  and  final  electron  momenta,      f      and  f     are  the  longitudinal 

i—i  t 

and  transverse  form  factors.      Regardless  of  the  model  used,    the  cross 
section  will  have  this  form  and  will  depend  on  these  two  form  factors. 
If  nuclear  recoil  is  taken  into  account,    the  cross-section  is   simply 

multiplied  by  a  factor 

1 


k      -  k    cos  9 


E  ' 

The  transverse  form  factor  can  be  easily  separated  out,    since  for 

o 
9    =    180     ,     the  coefficient  of  the  longitudinal  form  factor  goes  to  zero. 


Figure  1   shows  the  experimental  results  for  the  square  of  the 

12 

transverse  form  factor  for  electron  scattering  from  C        as  a  function 

of  momentum  transfer,    along  with  the  predictions  of  two  collective 
models  (the  Goldhaber  -  Teller  and  Steinwedel  -  Jensen  models)  and  two 
particle-hole  shell  models  (the  Brown  and  Gillet  -  Melkanoff  models). 
Figure  2  shows  the  experimental  results  for  the  transverse  form  factor 

for  0         and  the  predictions  of  the  Brown  and  Goldhaber  -  Teller  models. 

1  ?  1  A 

For  both  C    '"*  and   0       ,     the  observed  behavior  is  more  closely  predicted 

by  the  shell  models  than  by  either  of  these  two  collective  models.      Lewis 
and  Walecka  assert  that  this  is  conclusive  proof  that  the  collective  mode 
description  of  these  excited  states  is  not  adequate,    and  that  the  particle- 
hole  shell  model  description  is  correct.      The  purpose  of  this  paper  is 
to  develop  a  collective  model  which  can  predict  the  observed  behavior 
of  the  form  factor. 

In  the  Goldhaber  -  Teller  model,    the  neutrons  of  the  excited  nucleus 
are  assumed  to  move  against  the  protons  as  a  unit,    thus  giving  the 
effect  of  an  harmonic  oscillator,    with  the  restoring  force  being  due  to 
changes  in  the  overlap  of  the  protons  and  the  neutrons  [2].      The 
interaction  Hamiltonian  is  assumed  to  be 

12,  1  2 

H  =    -  P     /  u    +    -    uwQ      , 

where  w  is  the  frequency  of  oscillation,    u  =  NZM  /  A    is  the  reduced 
mass  of  the  system,     Q   is  the  coordinate  representing  the  separation 
of  the  center  of  mass  of  the  proton  system  from  the  center  of  mass  of 


the  neutron  system,    and  P    is  the  canonical  momentum  associated  with 
the  coordinate  Q.      P  and  Q  are  quantized  and  the  displacements  are 
assumed  to  be  small.      The  charge  and  current  densities  are  found  from 
the  motion  of  the  protons. 

In  the  shell  model  theories,    a  single  nucleon  is  assumed  to  be 
excited  from  a  closed  shell  by  the  interaction  with  the  incoming  electron. 
This  nucleon  then  interacts  with  the  resulting  "hole"  left  in  the   shell  it 
previously  occupied.      The  nucleon  and  the  hole  are  treated  as  point 
particles.     The  excited  state  of  the  nucleus  is  assumed  to  be  a  linear 
combination  of  the  single-particle  excitations,    and  the  charge,    current, 
and  magnetization  matrix  elements  between  the  single -particle  states 
are  found. 

In  this  paper,    a  collective  mode  is  assumed  and  a  semiclassical 
treatment  is  used  to  find  the  matrix  elements  of  the  Hamiltonian  de- 
scribing the  interaction  between  the  electron  and  the  nucleus.      The 
cross-section  is  then  found  to  depend  on  the  two  expected  form  factors. 
These  form  factors  are  then  calculated  for  a  nucleus  of  constant  density 
and  for  a  nucleus  whose  charge  distribution  is  the  Fermi  distribution. 
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2.       The   Model. 

'jollective  mode  assumed  is  an  harmonic  oscillation  of  the 
protons  and  neutrons  in  the  8 -direction  about  some  axis,    so  that  the 
char^<    density  of  the  excited  part  of  the  nucleus  is: 

p(r)     =     Z'e    f(r)     (1  +  a    cos  9    e  ) 

wh<  re   Z     is  the  number  of  protons  involved  in  the  excitation,    and 
Z    e    f(r)    is  the  ground  state  charge  density  of  this  portion  of  the  nucleus. 
It  is  assumed  that  the  neutrons1  motion  is  equal  and  opposite  to  that  of 
the  protons,    so  that  the  net  angular  momentum  is  zero.      At  this  point, 
no  specification  is  made  as  to  whether  the  entire  nucleus  or  a  portion 
only  is  excited,    although  it  may  later  be  shown  that  in  order  to  get  the 
desired  results  for  the  transverse  form  factor,    it  is  necessary  to  assume 
that  only  an  outer  core  is  excited,    with  an  inner  core  of  nucleons  re- 
maining inert,      uo  ,     the  frequency  of  the  oscillation,    is  equal  to  the 
energy  transferred  to  the  nucleus  by  the  interacting  electron.      The 
system  of  units  employed  in  this  paper  is  that  in  which    c  =  "h    -  1  . 
From  charge  conservation, 

*•''■-%    • 

The  assumed  mode  has 


J  =  Jeee      ' 

where     e      is  the  unit  vector  in  the   9    direction.      Thus,    the  charge 
8 
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conservation  equation  is; 


r   sine        d9 
This  leads  to: 


(J    sin0)     =     -iujaZ  e     f(r)      cose     e 


—     iaaaZ  e     r     f(r)      sine     e 


J(r)    =        2 
The  calculation  of  the  matrix  elements  H      is  done  semi-classically. 


The  matrix  element  is  [5] 


Hfi 


3 
J   (x)A   (x)     d  x   . 


It  is  assumed  that  the  nucleus  is  scattered  by  the  electric  field  of  the 
electron,    which  can  be  written  as 


4  Ilie  i  q    •   x    - 

— r e     -*       -    u 


A   (x)     =   — r e     -       —    u     y    u.    , 

*-._  I         LL      1 


q    V 

where  q  is  the  four -momentum  transfer  and  u„  and  u.  are  the  final 
—  f  i 

and  initial  plane-wave  electron  states  [5].      V  is  the  normalization 
volume  for  these  states.      It  is  further  assumed  that  J    (x)  for  the 
nucleus  can  be  written  as  just  the  classical  charge  and  current  densities 
for  the  model.      Thus: 


Hfi 


-r    ,    *   a     /    s       ,3              4  [lie         i  q  •  x         ,    .  -  3 

J    (x)  A   (x)     d    x    =   — x e   —     —   J   (x)  u    v    u.    d    x 


But    J   (x)    is  just   j    (x)    e  ,     and 


iq*x  iq*x         -icut 

e    —        —    =     e  e 
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so  what  must  be  calculated  is  the  Fourier  transform  of  the  current 

density: 

j    (q)    =      %iq         '   j    (x)    d  x     . 


1  I 


3.       Calculation  of  Fourier  Transforms. 

There  will  be  three  normal  modes  for  this  model:     one  for  which 

the  axis  of  vibration  is  parallel  to   q   (the  longitudinal  mode)  and  two  for 

which  the  axes  of  vibration  are  perpendicular  to    q   (the  transverse  modes), 

For  the  longitudinal  mode,    let  the  z-axis  denote  the  axis  of  vibration. 

Then: 

iqrcos  9 


J    (q)    =     C        r       f(r)     sin   9     cos9     cos0     e 


drd9d0 


=     0 


Likewise: 


p      3  2 

J     (q)     =     C        r       f(r)     sin   9     cos9     sin0      e 

y  j 


iqrcos  9 


drd9d0 


=     0 


Thus,    for  the  longitudinal  mode,    both  transverse  terms  are  zero. 
However: 


jz(q)     =      -C 


3    „.    ,        .3  iqrcos 9 

r      f(r)     sin  9      e  drd9d0 


snc 


f(r)     [  sinqr     -    qr     cos  (qr)  ]  dr 


=     4IliaZ/ea)      f_   (q) 
L-i 

where 

fL(q)     =    — J       I    [sinqr     -     qr     cos  (qr)  ]  dr    . 


The  fourth  component  of  j  (q)    is; 


JpM 


i  q  .  x       3 

e    —     —    d   x    =     ia 


,.    .       2  iqrcos9 

f(r)     r       sinfi     cos  9    e  drd9d0 


4  naZ  eq    f     (q)      . 
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Thus,    for  the  longitudinal  mode, 


j     (q)     =     4niaZ/eou     f     (q)     (0,   0,    1,    ^)    . 


For  the  transvrsc-  modes,    the  calculations  are  done  in  the  coordi- 
nat<       v  tt<  in    in  which  the  z-axis  is  the  axis  of  vibration  and  the  y-axis 
is   in  the  direction  of   q.      For  one  of  these  modes: 


i  p  (q) 


r „ .    ,     2  iqr      sine      sin0 

f  (r)  r         sine     cose      e  drded0    =     0 


jx(q)  C 


. ,    ,     3         .2  iqr      sine     sin0 

f  (r)  r        sin   Q  cose     cos0     e  drd8d0 


Jy(q] 


C 

0 


,,    ,     3         .    2  iqr      sine     sin0 

f(r)r        sin   Q  cos  Q     sin0      e  drded0 


j   (q) 

z 


3  3  iqr      sine      sin0 


C    j   f  (  r)  r  '      sin^G     c 

4n  c 


drded0 


f  (r)  [  qr    -   sin(qr)  ]  dr    =    -   2  flia  uu     f   (q) 


whe  r  e 


ft  ^q*   =  ~T    I  f  ^*  ^ sin  ^r)  -  qr ]  dr 


For  the  modes  in  which  the  axis  of  vibration  is  perpendicular  to    q  , 
the  only  nonzero  Fourier  transforms  are  those  in  the  directions  of 
the  axes  of  vibration.      The  overall  current  density  Fourier  transform 
can  be  written: 
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£(,)     =    4IIU..Z'.    [-    I4pl     ft(q),     -±6p2     f^q).    6p3    f^q). 

iS  6o*    fT  (q>        =     4nia«u  Z'e     fP(q) 

u)     P3      L  J  p, 

where  P  is  the  polarization  axis  for  the  excited  nucleus  (1,    2,    or  3, 
with  the  3 -axis  being  the  axis  parallel  to    q  )    . 
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4.        Calculation  of  Cross-Section. 

Th<    i  it'  raction  Hamiltonian  is  now: 


P  4nie       .P  - 

Hr:  5 J        (q)       Ur   Y      U; 


q    V 


2  2 

1 6  II      aZ'eoj  r  P  ,    .    - 

— f      (q)    u    y    u 

£--  p,  t     u-     i 


q    V 

The  cross-section  depends  on  the  square  of  the  Hamiltonian, 
averaged  over  initial  electron  spin  states  and  summed  over  final 
electron   spins  and  nuclear  polarizations.      This   is: 

2      „/    2         v  2  2 


t,iZ  l        r  /    16TI      aZ    e     uu      \        ..P.       - 

i,  f,  P  q    V 


I  i  ■     Mimmations  lead  to  a  trace  calculation,    as  follows: 

p  _  $  P  *  x 

£  (    i      (q)    u    v    u.  u.     v      v      u    f  )   = 

Tr     (    u.  v    u.  u.      y    y„  u,    ) 
LV    i   V     i     i       v  '4     fy 


P  *     P 
f  f 

p,ii  t     i  y     v         p, 


Let  us  now  examine  the  quantity  in  parentheses.      If  the    y       and 
y      were   reversed,    this  would  just  be 

Tr  (  u    v    u-  u-  "Y     ur    )  > 
which  is  [  5  ]  : 
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EE 


~    [k  k'     +    k  k  '    -    (k    *k'+m2)6  ]  . 


v  v 


From  the  commutation  relations  for  the  gamma  matrices: 


Tr    (5f^u.u.     YvY4uf)     = 


E1E£ 


[  k    k  '  *      +    k  *  k  '    -    (  k    *  k      +    m2)     (6         -26    ,6     „)] 


4      v4 


The  form  factor  matrix  is: 


F      =  t  f2  (q) 

vp,         4     t     VM/ 


1 

0 

0 

0 

0 

i 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

+   fL  (q) 


0    0         0 


0    0         0 


0    0         1  iq/uj 


2     2 

0    0    -iq/uu      q    /u) 


\   f 2  (q)    A'         +       f 2    (q)    A^ 

4     t  vjj,  L,  vp, 


There  are  two  traces  which  must  be  evaluated.      These  are; 


The  transverse  coefficient: 


Tt{B^aIJ=  4   sin2(e/2)  [i    --^r-    cos2(e/2) 


The  longitudinal  coefficient: 

2q2 
Tr  (Bv       ^J      =     "V      C°s2  (0/2) 

q 

where    e  is  the  angle  between  the  incoming  electron  momentum  k  and 
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the  outgoing  electron  momentum  k'.      The  square  of  the  Hamiltonian, 
averaged  and  summed  over  states,    is: 

l„,2  f  I6n2  aZ'e    N2     2  1  ,  ,        .    2 

lHi      (.t: — )  ™  iTkP  {2kk  sin  (9/2> 

T  V 


I    ,  2kk  2    ,     . _,   ~l    ,2    .    , 

X         1      +    —r      cos       (6/2)    J    ft     (q) 


A,    l      '         4 

,      4kk     q  2  Z  1 

-TTT-     cos     <6/2)    fL  (q)  }   • 

q     a) 


The  amplitude    a    of  the  oscillation  can  be  calculated  (classically) 
from  energy  considerations.      At  the  time  during  the  oscillation  such 
that  the  charge  density  of  the  nucleus  is  its  ground-state  configuration 
(  uu  t  (n  +  —  )  II  )  ,     all  the  energy  transferred  to  the  nucleus  is  in  the 

kinetic  energy  of  the  neutrons  and  protons.      The  velocity  distribution 
of  the  excited  protons  is 

v   (r)     =    "j  (r)/  Z'ef(r)    . 

This  ratio  is    v(r)     =    —  oiar      sinG  .      The  mass  density  of  the  excited 
protons  is    Z    M    f  (r).      Thus,    the  kinetic  energy  distribution  of  the 
excited  protons  is: 

^   MZ'   f(r)     v2(r)     =    i   Z'M    f(r)    uu       a    r        sin  6    . 
2  o 

Since  half  of  the  transferred  energy  will  be  in  the  protons  (with  the 

other  half  being  in  the  kinetic  energy  of  the  neutrons): 
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,    2    2 

u)  MZ   ou    a         ^ 


2  2         2 

f(r)     r       sine     r       sin  9      d9    dr      d0 


This  gives 


211  Z'MooI 


where    I    = 


Jf(r) 


4 

r       dr  . 


The  cross -section  is,    from  the  well-known  Golden  Rule  of  Fermi, 


da     -       L  211 

final 
states 


Hi2    5    (E  -E.)     -i- 
f l  '  f        l         flux 


The  flux  is  just  k/  VE.  .      The  summation  over  final  states  is  just 

3         ,  r 

an  integral  over    d    k     for  the  case  at  hand.      The    o-function  is  just 
6    (k  -  k'-  w  )  .      Thus: 

da     =    J- 


V         VEi 


(2n) 


H     |         £  (k  -  k'  -  cu)    k'  '    dk'  dQ 


2  ,2 

V   E.  k' 

2 I  Hr  I        dQ 

(2n)   k  " 


k  -  k      =    u) 


Finally,    putting    Ih|       into  this  equation,    we  obtain: 


2    , 

da  a    k 

du>  4 

q   k  2 


-^ I    2kk'     sin2(Q/2)     f      [l+-^-     cos2(Q/2) 

n  z'mi     l  t  L       q2  J 


4kk    q 

~~ 2      2~ 

q    uj 


fL       cos     (G/2)   |       . 
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Note  that  these  form  factors  must  be  multiplied  by  a  constant 


ID 


2n2  z'mi 


in  order  to  compare  them  with  the  form  factors  shown  in  Figures   1 
and  2. 
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5.       Uniform  Shell  Model. 

If  the  nuclear  charge  density  is  assumed  to  be  uniform,  with  the 
nucleons  within  a  shell  of  inner  radius  r  and  outer  radius  r  being 
excited  by  the  interaction,    we  have: 

fL(q)     =    —      Tsin(qr2)    -    qr2      cosfqr^)     -    [sirifqi^)     -    q^      cosfqr^jj 

q 

ft    (q)     =     —     rsin(qr2)    -    qr2    -    [sintqi^)    -    q^] 

q 

Figure   3  shows  the  results  for  the  square  of  the  transverse  form 
factor  versus  momentum  transfer  for  0       ,     with  the  assumption  that  the 

entire  nucleus  is  excited   (r      -  0)  ,     and  the  same  form  factor  with  the 

4 
assumption  that  an  inner  core  consisting  of  a   He      nucleus  remains 

inert.      The  figures  used  for  the  outer  radii  for  these  calculations  were 

3.  41  fermi  for  0         and  2.  08  fermi  for  He        [3]. 

If  one  assumes  an  even  thinner  shell  is  excited,    the  form  factor  is 

reduced  by  orders  of  magnitude  and  falls  off  to  zero  with  increasing    q 

even  faster.      Figure  4  shows  the  results  for  a  shell  of  thickness  0.  1 

fermi.      No  matter  how  thin  the  shell,    the  uniform  shell  model  cannot 

give  the  observed  form  factor  behavior  with    q   for  a  reasonable  outer 

radius. 
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6.         I  h<     :  ■  i  mi  Distribution  Model. 

Now  let  us  assume  that  the  charge  density  function  is  just  the 
Fermi  distr  ibut  ion: 
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integral  can  be  done  by  breaking  it  up  into  two  parts  and  writing 
ea<  h  pari    as  an  infinite'   .scries,    as  follows: 
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The    result   for  the  form  factor  is: 
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For  excitation  of  a  shell,    the  form  factor  is 

ft(q)   =   ita  (q)    -   £  (q) 

where    a    is  the  outer  radius  of  the  shell  and  b    is  the  inner  radius. 

The  square  of  the  transverse  form  factor  for  0         is  shown  in  Figure   5, 

using   a  =  3.  64    and  these  combinations  for  b  and  a:    b  =  3.  0  , 

a    =    1;    b  =  3.6,     a  =  1  ;     and  b  =  3.  0  ,     ot  =    20  .      Again,    the  thinner 

the  shell,    the  smaller  the  form  factor  and  the  faster  it  drops  off  with 

increasing   q.      For  large  a,    the  square  of  the  form  factor  is  smaller 

and  does  not  drop  off  as  fast.      In  fact,    for  very  large    a,     one  gets  just 

2 

the  uniform  shell.      Note  that  the  f        curve  for    a  -  20    actually  crosses 

the  oi   -  1    curve.      This  leads  one  to  suspect  that  if  one  assumes  an  outer 
boundary  for  the  excited  shell  with,    say,     a   -   1    and  an  inner  boundary 

with  very  large    a  ,    for  a  thin  enough  shell  one  could  get  a  curve  for 

.         2 
f  which  has  the  desired  minimum  at  about    q  -  0.  3  . 

Figure  6  shows  the  results  for  three  trials  of  this  composite  model. 

For  a  thick  shell    (a  =  3.  64    and   b  -  2.  0)    with   a  -   1  ,     the  form  factor  is 

larger  than  for  the  same  shell  with   a  =  10.      The  form  factor  becomes 

larger  for  a  thinner  shell.      However,    for  the  combinations  of  a,    b,    and  a 

used  here,    the  form  factor  does  not  get  the  desired  behavior. 
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7.        Conclusions. 

t     None  of  the  assumed  charge  distributions  gives  the  desired  behavior 

for  the  form  factor.      This  model  gives  somewhat  the  same  results  as 

i  i  2 

the  Goldhaber  -Teller  model;  that  is,    the  function       f  (q)  is  a 

motonically  decreasing  function  of  q   within  the  range  of  interest. 
Further  work  is  necessary  on  the  composite  model  before  this  model 
can  be  completely  discarded,    however.      For  instance,    if  one  assumes 
the  nucleus  as  a  whole  has  a  Fermi  charge  distribution  and  that  the  inert 
core  is  a  uniform  shell  with  a  radius  equal  to  the  radius  of  the  half- 
height  of  the  Fermi  distribution,    what  one  has  is  essentially  a  classical 
treatment  of  the  particle-hole  theory.      Figure  7  illustrates  the 
"particle"  and  the  "hole"  associated  with  this  model.      This  idea  should 
be  pursued  further.      There  is  a  strong  possibility  that  it  will  give  re- 
sults similar  to  those  of  the  successful  Brown  particle-hole  theory. 
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